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Abstract We discuss the different physical processes that are important to under- 
stand the thermal X-ray emission and absorption spectra of the diffuse gas in clusters 
of galaxies and the warm-hot intergalactic medium. The ionisation balance, line and 
continuum emission and absorption properties are reviewed and several practical ex- 
amples are given that illustrate the most important diagnostic features in the X-ray 
spectra. 
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1 Introduction 



Thermal X-ray radiation is an important diagnostic tool for studying cosmic sources 
where high-energy processes are important. Examples are the hot corona of the Sun 
and of stars, Solar and stellar flares, supernova remnants, cataclysmic variables, ac- 
cretion disks in binary stars and around black holes (Galactic and extragalactic), the 
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diffuse interstellar medium of our Galaxy or external galaxies, the outer parts of Active 
Galactic Nuclei (AGN), the hot intracluster medium, the diffuse intercluster medium. 
In all these cases there is thermal X-ray emission or absorption. 

In the present paper we focus upon the properties of the X-ray spectrum in the 
diffuse gas within and between galaxies, clusters and the large scale structures of the 
Universe. As this gas has very low densities, the level populations in the atoms are not 
governed by Saha-like equations, but instead most of the atoms will be in or near the 
ground state. This simplifies the problem considerably. Furthermore, in these tenuous 
media we need to consider only gas with small or moderate optical depth; the radiative 
transport is therefore simple. Stellar coronae have similar physical conditions, but there 
the densities are higher so that in several cases the emergent spectrum has density- 
dependent features. Due to the low densities in our sources, these density effects can 
be ignored in most cases; however, for the lowest density gas photoionisation effects 
must be taken into account. 

We show in this paper that it is possible to derive many different physical parame- 
ters from an X-ray spectrum: temperature, density, chemical abundances, plasma age, 
degree of ionisation, irradiating continuum, geometry etc. 

The outline of this paper is as follows. First, we give a brief overview of atomic 
structure (Sect. [2]). We then discuss a few basic processes that play an important role 
for the thermal plasmas considered here (Sect. |3}. For the proper calculation of an 
X-ray spectrum, three different steps should be considered: 

1. the determination of the ionisation balance (Sect.|4)|, 

2. the determination of the emitted spectrum (Sect. [5}, 

3. possible absorption of the photons on their way towards Earth (Sect. [6]). 

We then briefly discuss issues like the Galactic foreground emission (Sect. [7} , plasma 
cooling (Sect. [SJ, the role of non-thermal electrons (Sect. [9]), and conclude with a section 
on plasma modelling (Sect. [TO)) . 

2 A short introduction to atomic structure 

2.1 The Bohr atom 

The electrons in an atom have orbits with discrete energy levels and quantum numbers. 
The principal quantum number n corresponds to the energy I n of the orbit (in the 
classical Bohr model for the hydrogen atom the energy I n = E-R_n~ 2 with En = 13.6 eV 
the Rydberg energy), and it takes discrete values n =1, 2, 3, . . .. An atomic shell consists 
of all electrons with the same value of n. 

The second quantum number £ corresponds to the angular momentum of the elec- 
tron, and takes discrete values £ < n. Orbits with £ =0, 1, 2, 3 are designated as s, p, 
d, and f orbits. A subshell consists of all electrons with the same value of n and £; they 
are usually designated as Is, 2s, 2p, etc. 

The spin quantum number s of an electron can take values s — ±1/2, and the 
combined total angular momentum j has a quantum number with values between 
£ — 1/2 (for £ > 0) and £+ 1/2. Subshells are subdivided according to their j-value and 
are designated as n£j. Example: n = 2, £ = 1, j = 3/2 is designated as 2p 3/ / 2 - 

There is also another notation that is commonly used in X-ray spectroscopy. Shells 
with n = 1, 2, 3, 4, 5, 6 and 7 are indicated with K, L, M, N, O, P, Q. A further 
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subdivision is made starting from low values of I up to higher values of t and from low 
values of j up to higher values of j: 

Is 2s 2p 1/2 2p 3/2 3s 3p 1/2 3p 3/2 3d 3/2 3d 5/2 4s etc. 

K L/ L77 L777 M/ M77 M777 M IV My N7 
222 422 4 4 62 

The third row in this table indicates the maximum number of electrons that can be 
contained in each subshell. This so-called statistical weight is simply 2j + 1. 

Atoms or ions with multiple electrons in most cases have their shells filled starting 
from low to high n and £. For example, neutral oxygen has 8 electrons, and the shells 
are filled like ls 2 2s 2 2p 4 , where the superscripts denote the number of electrons in each 
shell. Ions or atoms with completely filled subshells (combining all allowed j-values) 
are particularly stable. Examples are the noble gases: neutral helium, neon and argon, 
and more general all ions with 2, 10 or 18 electrons. In chemistry, it is common practice 
to designate ions by the number of electrons that they have lost, like 2+ for doubly 
ionised oxygen. In astronomical spectroscopic practice, more often one starts to count 
with the neutral atom, such that 2+ is designated as O in. As the atomic structure 
and the possible transitions of an ion depend primarily on the number of electrons it 
has, there are all kinds of scaling relationships along so-called iso-electronic sequences. 
All ions on an iso-electronic sequence have the same number of electrons; they differ 
only by the nuclear charge Z. Along such a sequence, energies, transitions or rates are 
often continuous functions of Z. A well known example is the famous Is— 2p triplet 
of lines in the helium iso-electronic sequence (2-electron systems), e.g. Cv, Nvi and 
Ovn. 



2.2 Level notation in multi-electron systems 

For ions or atoms with more than one electron, the atomic structure is determined 
by the combined quantum numbers of the electrons. These quantum numbers have 
to be added according to the rules of quantum me chanics. We will not go into detail 
here, but refer to textbooks such as lHerzbersI ( 19441 ). The determination of the allowed 



quantum states and the transitions between these states can be rather complicated. 
Important here is to know that for each electron configuration (for example Is 2p) there 
is a number of allowed terms or multiplets, designated as 2S+1 L, with 5 the combined 
electron spin (derived from the individual s values of the electrons), and L represents 
the combined angular momentum (from the individual £ values). For L one usually 
substitutes the alphabetic designations similar to those for single electrons, namely S 
for L — 0, P for L = 1, etc. The quantity 25 + 1 represents the multiplicity of the 
term, and gives the number of distinct energy levels of the term. The energy levels of 
each term can be distinguished by J, the combined total angular momentum quantum 
number j of the electrons. Terms with 25*+ 1 equals 1, 2 or 3 are designated as singlets, 
doublets and triplets, etc. 

For example, a configuration with two equivalent p electrons (e.g., 3p 2 ), has three 
allowed terms, namely 1 S, *D and 3 P. The triplet term 3 P has 25 + 1 = 3 hence 5=1 
and L — 1 (corresponding to P) , and the energy levels of this triplet are designated as 
3 Pq, Pi and 3 P 2 according to their J values 0, 1 and 2. 
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Fig. 1 Energy levels of atomic subshells for neutral atoms. 



Table 1 Binding energies E, corresponding wavelengths A and photoionisation cross sections 
<r at the edges of neutral atoms (nuclear charge Z) for abundant elements. Only edges in the 
X-ray band are given. Note that 1 Mbarn = 10~ 22 m 2 . 
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(Mbarn) 


Is shell: 








2s shell: 








H 


1 


13.6 


911.8 


6.29 


O 


8 


16.6 


747.3 


1.37 


He 


2 


24.6 


504.3 


7.58 


Si 


14 


154 


80.51 


0.48 


C 


6 


288 


43.05 


0.96 


s 


16 


232 


53.44 


0.37 


N 


7 


403 


30.77 


0.67 


Ar 


18 


327 


37.97 


0.29 


O 


8 


538 


23.05 


0.50 


Ca 


20 


441 


28.11 


0.25 


Ne 


10 


870 


14.25 


0.29 


Fe 


26 


851 


14.57 


0.15 


Mg 


12 


1308 


9.48 


0.20 


Ni 


28 


1015 


12.22 


0.13 


Si 


14 


1844 


6.72 


0.14 
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16 


2476 


5.01 


0.096 
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16 
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Ar 


18 
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3.87 
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Ar 


18 


251 
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1.45 


Ca 


20 
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3.07 


0.060 


Ca 


20 


353 


35.16 


0.86 


Fe 


26 


7117 


1.74 


0.034 


Fe 


26 


726 


17.08 


0.42 


Ni 


28 


8338 


1.49 


0.029 


Ni 


28 


877 


14.13 


0.35 
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18 
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20 
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35.53 
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Fe 


26 


713 


17.39 


0.86 












Ni 


28 


860 


14.42 


0.73 
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2.3 Binding energies 

The binding energy / of K-shell electrons in neutral atoms increases approximately as 
I ~ Z 1 with Z being the nuclear charge (see Table [1] and Fig. [1} . Also for other shells 
the energy increases strongly with increasing nuclear charge Z. 



Table 2 Binding energies E, corresponding wavelengths A and photoionisation cross sections 
<t at the K-edges of oxygen ions. 



ion 


E 


A 


a 




(eV) 


(A) 


(Mbarn) 
(10^ 22 m 2 ) 


Oi 


544 


22.77 


0.50 


On 


565 


21.94 


0.45 


Om 


592 


20.94 


0.41 


Oiv 


618 


20.06 


0.38 


Ov 


645 


19.22 


0.35 


Ovi 


671 


18.48 


0.32 


Ovn 


739 


16.77 


0.24 


O viii 


871 


14.23 


0.10 



For ions of a given element the ionisation energies decrease with decreasing ionisa- 
tion stage: for lowly ionised ions, a part of the Coulomb force exerted on an electron 
by the positively charged nucleus is compensated by the other electrons in the ion, 
thereby allowing a wider orbit with less energy. An example is given in Table [5] 



2.4 Abundances 

With high spectral resolution and sensitivity, optical spectra of stars sometimes show 
spectral features from almost all elements of the Periodic Table, but in practice only 
a few of the most abundant elements show up in X-ray spectra of cosmic plasmas. 
In several situations the abundances of the chemical elements in an X-ray source are 
similar to (but not necessarily equal to) the abundances for the Sun or a constant 
fraction of that. There have been many significant changes over the last several years 
in the ado pted set of standard cosmi c abundances. A still often used set of abundances 
is that of An ders fc Grevessel l|l989h . but a more recent one is the set of proto-solar 



abundances of iLoddersl (|2003r ). that we list in Table [3] for a few of the key elements. 

In general, for absorption studies the strength of the lines is mainly determined 
by atomic parameters that do not vary much along an iso-electronic sequence, and 
the abundance of the element. Therefore, in the X-ray band the oxygen lines are the 
strongest absorption lines, as oxygen is the most abundant metal. The emissivity of 
ions in the X-ray band often increases with a strong power of the nuclear charge. For 
that reason, in many X-ray plasmas the strongest iron lines are often of similar strength 
to the strongest oxygen lines, despite the fact that the cosmic abundance of iron is only 
6 % of the cosmic oxygen abundance. 
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Table 3 Proto-solar abundances for the 15 m ost common ch emical elements. Abundances A 
are given with respect to hydrogen. Data from Loddcrs (2003). 

Element abundance Element abundance Element abundance 



H 


= 1 


Ne 


89.1 x 10" 6 


S 


18.2 x 10~ 6 


He 


0.0954 


Na 


2.34 x 10~ 6 


Ar 


4.17 x 10~ 6 


C 


288 x 10" 6 


Mg 


41.7 x 10~ 6 


Ca 


2.57 x 10~ 6 


N 


79.4 x 10~ 6 


Al 


3.47 x 10~ 6 


Fe 


34.7 x 10~ 6 





575 x 10" 6 


Si 


40.7 x 10" 6 


Ni 


1.95 x 10~ 6 



3 Basic processes 

3.1 Excitation processes 

A bound electron in an ion can be brought into a higher, excited energy level through a 
collision with a free electron or by absorption of a photon. The latter will be discussed 
in more detail in Sect. [6] Here we focus upon excitation by electrons. 

The cross section Qij for excitation from level i to level j for this process can be 
conveniently parametrised by 

where U = E^/E with E^ the excitation energy from level i to j, E the energy of the 
exciting electron, the Rydberg energy (13.6 eV), ap the Bohr radius and ittj the 
statistical weight of the lower level i. The dimensionless quantity Q(U) is the so-called 
collision strength. For a given transition on an iso-electronic sequence, fl(U) is not a 
str ong function of the atomic number Z, or may be even almost independent of Z. 

iMewel (|l972l ) introduced a convenient formula that can be used to describe most 
collision strengths, written here as follows: 

(2(U)=A+^ + £L + ™+FlnU, (2) 

where A, B, C, D and F are parameters that differ for each transition. The expression 
can be integrated analytically over a Maxwellian electron distribution, and the result 
can be expressed in terms of exponential integrals. We only mention here that the total 
excitation rate Sij (in units of m^ 3 s^ 1 ) is given by 

^8.62X1CT 12 fiW- 1 ' 2 *-* , ( 3 ) 

Wi 

with y = Eij/kT and Q{y) is the Maxwellian-averaged collision strength. For low 
temperatures, y > 1 and D(y) = A + B + C + 2D, leading to S i3 ~ T~ 1 / 2 e" y . 
The excitation rate drops exponentially due to the lack of electrons with sufficient 
energy. In the other limit of high temperature, y < 1 and fi(y) = — Flny and hence 
~T- 1 /2]ni/. 

Not all transitions have this asymptotic behaviour, however. For instance, so-called 
forbidden transitions have F = and hence have much lower excitation rates at high 
energy. So-called spin- forbidden transitions even have A = B = F = 0. 
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In most cases, the excited state is stable and the ion will decay back to the ground 
level by means of a radiative transition, either directly or through one or more steps 
via intermediate energy levels. Only in cases of high density or high radiation fields, 
collisional excitation or further radiative excitation to a higher level may become im- 
portant, but for typical cluster and ISM conditions these processes are not important 
in most cases. Therefore, the excitation rate immediately gives the total emission line 
power. 



3.2 Ionisation processes 
3.2.1 Collisional ionisation 

Collisional ionisation occurs when during the interaction of a free electron with an atom 
or ion the free electron transfers a part of its energy to one of the bound electrons, which 
is then able to escape from the ion. A necessary condition is that the kinetic energy E 
of the free electron must be larger than the binding energy / of the atomic shell from 
which the bound electron escapes. A formula that gives a correct order of magnitude 
estimate of the cross section a of this process and that has t he proper a symptotic 
behaviour (first calculated by Bethe and Born) is the formula of lLotzl ( 19681 ): 



_ an s lnQg/J) , . 

EI ' 1 ' 

where n s is the number of electrons in the shell and the normalisation a — 4.5 x 
10 -24 m 2 keV 2 . This equation shows that high-energy electrons have less ionising power 
than low-energy electrons. Also, the cross section at the threshold E = I is zero. 

The above cross section has to be averaged over the electron distribution (Maxwellian 
for a thermal plasma). For simple formulae for the cross section such as Q the integra- 
tion can be done analytically and the result can be expressed in terms of exponential 
integrals. We give here only the asymptotic results for Cdi, the total number of direct 
ionisations per unit volume per unit time: 

kT«I: Cd] ,(^)W|^, (5) 
V y/irmZ / I z 

and 

kT»/ : c ^(?^y^Mw/i) 

\ ^nT c ) J^/kT 

For low temperatures, the ionisation rate therefore goes exponentially to zero. This 
can be understood simply, because for low temperatures only the electrons from the 
exponential tail of the Maxwell distribution have sufficient energy to ionise the atom 
or ion. For higher temperatures the ionisation rate also approaches zero, but this time 
because the cross section at high energies is small. 

For each ion the direct ionisation rate per atomic shell can now be determined. The 
total rate follows immediately by adding the contributions from the different shells. 
Usually only the outermost two or three shells are important. That is because of the 
scaling with I~ 2 and 7" 1 in ([5} and ([6]), respectively. 
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Energy (keV) Energy (keV) 

Fig. 2 Photoionisation cross section in barn (10 -28 m — 2 ) for Fel (left) and Fexvi (right). 
The p and d states (dashed and dotted) have two lines each because of splitting into two 
sublevels with different j (see Sect, |27TT| . 



3.2.2 Photoionisation 



This process is very similar to collisional ionisation. The difference is that in the case 
of photoionisation a photon instead of an electron is causing the ionisation. Further, 
the effective cross section differs from that for collisional ionisation. As an example 
Fig. [5] shows the cross section for neutral iron and Na-like iron. For the more highly 
ionised iron the so-called "edges" (corresponding to the ionisation potentials 7) occur at 
higher energies than for neutral iron. Contrary to the case of collisional ionisation, the 
cross section at threshold for photoionisation is not zero. The effective cross section 
just above the edges sometimes changes rapidly (the so-called Cooper minima and 
maxima) . 

Contrary to collisional ionisation, all the inner shells now have the largest cross 
section. For the K-shell one can approximate for E > I 

a(E) ~ a (I/E) 3 . (7) 



For a given ionising spectrum F(E) (photons per unit volume per unit energy) the 
total number of photoionisations follows as 



Cpi = cJ n- l a(E)F{E)AE. (8) 


For hydrogenlike ions one can write: 

where n is the principal quantum number, a the fine structure constant and ao the Bohr 
radius. The Gaunt factor g(E,n) is of order unit y and varies only s lowly as a function 
of E. It has been calculated and tabulated by iKarzas fc Latter] 1 196lh . The above 



equation is also applicable to excited states of the atom, and is a good approximation 
for all excited atoms or ions where n is larger than the corresponding value for the 
valence electron. 
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3.2.3 Compton ionisation 



Scattering of a photon on an electron generally leads to energy transfer from one of 
the particles to the other. In most cases only scattering on free electrons is considered. 
But Compton scattering also can occur on bound electrons. If the energy transfer from 
the photon to the electron is large enough, the ionisation potential can be overcome 
leading to ionisation. This is the Compton ionisation process. 
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1 

Photon energy (keV) 
Fig. 3 Compton ionisation cross section for neutral atoms of H, He, C, N, O, Ne and Fc. 



In the Thomson limit the differential cross section for Compton scattering is given 



by 



dcr 3<T T 2 n-, 

d72 = l6^ (1 + COS 6) > 



(10) 



with 6 the scattering angle and <tt the Thomson cross section (6.65 x 10 m " 
The energy transfer AE during the scattering is given by (E is the photon energy): 



AE = 



E 2 {1 - cosfl) 
m c c 2 + E(l - cos9) ' 



(11) 



Only those scatterings where AE > I contribute to the ionisation. This defines a 
critical angle 9 C , given by: 

cosg c = l- „^ c2 /„ . (12) 



E 2 - IE' 

For E ^> I we have o(E) — > or (all scatterings lead in that case to ionisation) and 
further for 6 C — > n we have o(E) — > 0. Because for most ions / <§C m e c 2 , this last 
condition occurs for E ~ \J Im e c 2 /2 I. See Fig. [3] for an example of some cross 



sections. In general, Compton ionisation is important if the ionising spectrum contains 
a significant hard X-ray contribution, for which the Compton cross section is larger 
than the steeply falling photoionisation cross section. 
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3.2.4 Autoionisation and fluorescence 

As we showed above, interaction of a photon or free electron with an atom or ion 
may lead to ionisation. In particular when an electron from one of the inner shells is 
removed, the resulting ion has a "vacancy" in its atomic structure and is unstable. Two 
different processes may occur to stabilise the ion again. 

The first process is fluorescence. Here one of the electrons from the outer shells 
makes a radiative transition in order to occupy the vacancy. The emitted photon has 
an energy corresponding to the energy difference between the initial and final discrete 
states. 

The other possibility to fill the gap is auto-ionisation through the Auger process. 
In this case, also one of the electrons from the outer shells fills the vacancy in the 
lower level. The released energy is not emitted as a photon, however, but transferred 
to another electron from the outer shells that is therefore able to escape from the ion. 
As a result, the initial ionisation may lead to a double ionisation. If the final electron 
configuration of the ion still has holes, more auto-ionisations or fluorescence may follow 
until the ion has stabilised. 







~~~~ K 




L3„.'.': 




/ Us' 







0.4 



50 



100 




2 4 6 

Number of electrons 



Fig. 4 Left panel: Fluorescence yield u as a function of atomic number Z for the K and 
L shells. Right panel: Distribution of number of electrons liberated after the initial removal 
of an electron from the K-s hell (solid line) or Lj s hell (dotted line), including the original 
photo-electron, for Fel; after Kaastra & Mcwei \ 19931) . 



In Fig. U] the fluorescence yield u (probability that a vacancy will be filled by 
a radiative transition) is shown for all elements. In general, the fluorescence yield 
increases strongly with increasing nuclear charge Z, and is higher for the innermost 
atomic shells. As a typical example, for Fel a K-shell vacancy has lu = 0.34, while an 
Lr-shell vacancy has ui = 0.001. For Ol these numbers are 0.009 and 0, respectively. 

3.2.5 Excitation- Autoionisation 

In Sect. 13.2.11 we showed how the collision of a free electron with an ion can lead to 
ionisation. In principle, if the free electron has insufficient energy (E < I), there will 
be no ionisation. However, even in that case it is sometimes still possible to ionise, but 
in a more complicated way. The process is as follows. The collision can bring one of 
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the electrons in the outermost shells in a higher quantum level (excitation). In most 
cases, the ion will return to its ground level by a radiative transition. But in some 
cases the excited state is unstable, and a radiationless Auger transition can occur (see 
Sect. [37274]) . The vacancy that is left behind by the excited electron is being filled by 
another electron from one of the outermost shells, while the excited electron is able to 
leave the ion (or a similar process with the role of both electrons reversed). 



Phase 1: Phase 2: Final state: 




Fig. 5 Left panel: the Excitation-Autoionisation process for a Li-like ion. Right panel: col- 
lisional ionisation cross section for Fexvi. The contribution from direc t collisional ionisation 
(DI) and excitation-autoionisation (EA) are indicated. Adapted from lArnaud fc R othcnflue 
jl985ft . 



Because of energy conservation, this process only occurs if the excited electron 
comes from one of the inner shells (the ionisation potential barrier must be taken 
anyhow). The process is in particular important for Li-like and Na-like ions, and for 
several other individual atoms and ions. As an example we treat here Li-like ions (see 
Fig. [5] left panel). In that case the most important contribution comes from a Is— 2p 
excitation. 



3.3 Recombination processes 
3.3.1 Radiative recombination 

Radiative recombination is the reverse process of photoionisation. A free electron is 
captured by an ion while emitting a photon. The released radiation is the so-called 
free-bound continuum emission. It is relatively easy to show that there is a simple 
relation between the photoionisation cross section a^f(E) and the recombination cross 
section crft,, namely the Milne-relation: 

E 2 g n a M (E) 

O-fb(w) = = 2~ 



where g n is the statistical weight of the quantum level into which the electron is 
captured (for an empty shell this is g n = 2n 2 ). By averaging over a Maxwell distribution 
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one gets the recombination-coefficient to level n: 

oo 

Rn = nin c J u/(t;)(j fb (u)dii. (14) 


Of course there is energy conservation, so E = ^m c v 2 + I. 

It can be shown that for the photoionisation cross section © and for g = 1, constant 
and g n = 2n 2 : 

3V3m c kTZ 2 kTm cC 3 W ; V ; 

With the asymptotic relations for the exponential integrals it can be shown that 

kT</: i?„~T-5, (16) 
kT > J : R n ~ ln(J/kT)r~ 3/2 . (17) 

Therefore for T — > the recombination coefficient approaches infinity: a cool plasma 
is hard to ionise. For T — > oo the recombination coefficient goes to zero, because of the 
Milne relation (v — > oo) and because of the sharp decrease of the photoionisation cross 
section for high energies. 

As a rough approximation we can use further that I ~ (Z/n) 2 . Substituting this 
we find that for kT <C I (recombining plasmas) R n ~ n _1 , while for kT 3> / (ionising 
plasmas) R n ~ n -3 . In recombining plasmas in particular many higher excited levels 
will be populated by the recombination, leading to significantly stronger line emission. 
On the other hand, in ionising plasmas (such as supernova remnants) recombination 
mainly occurs to the lowest levels. Note that for recombination to the ground level 
the approximation (|15[1 cannot be used (the hydrogen limit), but instead one should 
use the exact photoionisation cross section of the valence electron. By adding over all 
values of n and applying an approximation Iseatonl 1 1959^ found for the total radiative 
recombination rate «rr (in units of m _3 s~ 1 ): 

a RR = ^R n = 5.197 x 10~ 20 ZA 1/2 {0.4288 + 0.5 In A + 0.469A _1/3 } (18) 

n 

with A = EftZ 2 /kT and En the Rydberg energy (13.6 eV). Note that this equation 
only holds for hydrogen-like ions. For other ions usually an analytical fit to numerical 
calculations is used: 

Q RR ~ T~ n (19) 
where the approximation is strictly speaking only valid for T near the equilibrium 
concentration. The approximations (|16p and (I17|) raise suspicion that for T — > or 
T — *• oo (|19p could be a poor choice. 

The captured electron does not always reach the ground level immediately. We 
have seen before that in particular for cool plasmas (kT I) the higher excited levels 
are frequently populated. In order to get to the ground state, one or more radiative 
transitions are required. Apart from cascade corrections from and to higher levels the 
recombination line radiation is essentially given by (|15[) . A comparison of recombina- 
tion with excitation tells that in particular for low temperatures (compared to the line 
energy) recombination radiation dominates, and for high temperatures excitation ra- 
diation dominates. This is also one of the main differences between photoionised and 
collisionally ionised plasmas, as photoionised plasmas in general have a low temperature 
compared to the typical ionisation potentials. 
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3.3.2 Dielectronic recombination 

This process is more or less the inverse of excitation-autoionisation. Now a free electron 
interacts with an ion, by which it is caught (quantum level n"l") but at the same time 
it excites an electron from (n£) — > (n £ ). The doubly excited state is in general not 
stable, and the ion will return to its original state by auto-ionisation. However there is 
also a possibility that one of the excited electrons (usually the electron that originally 
belonged to the ion) falls back by a radiative transition to the ground level, creating 
therefore a stable, albeit excited state (n"i") of the ion. In particular excitations with 
I = £ + 1 contribute much to this process. In order to calculate this process, one should 
take account of many combinations (n 1 £'){n" £"). 

The final transition probability is often approximated by 

«DR = ^- T °/ T (l + Be~ T ^ T ) (20) 

where A, B, Tq and T\ are adjustable parameters. Note that for T — > oo the asymp- 
totic behaviour is identical to the case of radiative recombination. For T — > however, 
dielectronic recombination can be neglected; this is because the free electron has in- 
sufficient energy to excite a bound electron. Dielectronic recombination is a dominant 
process in the Solar corona, and also in other situations it is often very important. 

Dielectronic recombination produces more than one line photon. Consider for ex- 
ample the dielectronic recombination of a He-like ion into a Li-like ion: 

c + Is 2 -> Is2p3s -> ls 2 3s + hvi -> ls 2 2p + hv 2 -> ls 2 2s + hu 3 (21) 

The first arrow corresponds to the electron capture, the second arrow to the stabil- 
ising radiative transition 2p^ls and the third arrow to the radiative transition 3s— >2p 
of the captured electron. This last transition would have also occurred if the free elec- 
tron was caught directly into the 3s shell by normal radiative recombination. Finally, 
the electron has to decay further to the ground level and this can go through the nor- 
mal transitions in a Li-like ion (fourth arrow). This single recombination thus produces 
three line photons. 

Because of the presence of the extra electron in the higher orbit, the energy hvi of 
the 2p— >ls transition is slightly different from the energy in a normal He- like ion. The 
stabilising transition is therefore also called a satellite line. Because there are many 
different possibilities for the orbit of the captured electron, one usually finds a forest 
of such satellite lines surrounding the normal 2p^ls excitation line in the He- like ion 
(or analogously for other iso-electronic sequences) . Fig. [S] gives an example of these 
satellite lines. 

3.4 Charge transfer processes 

In most cases ionisation or recombination in collisionally ionised plasmas is caused by 
interactions of an ion with a free electron. At low temperatures (typically below 10 5 K) 
also charge transfer reactions become important. During the interaction of two ions, 
an electron may be transferred from one ion to the other; it is usually captured in an 
excited state, and the excited ion is stabilised by one or more radiative transitions. 
As hydrogen and helium outnumber by at least a factor of 100 any other element (see 
Table [3} , in practice only interactions between those elements and heavier nuclei are 
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6.6 6.65 
Energy (keV) 

Fig. 6 Spectrum of a plasma in collisional ionisation equilibrium with kT = 2 keV, near the 
Fe-K complex. Lines are labelled using the most common designations in this field. The Fe XXV 
"triplet" consists of the resonance line (w), intercombination line (actually split into x and y) 
and the forbidden line (z). All other lines are satellite lines. The labelled satellites are lines 
from Fexxiv, most of the lines with energy below the forbidden (z) line are from Fexxm. 
The relative intensity of these satellites is a strong indicator for the physical conditions in the 
source. 



important. Reactions with Hi and He I lead to recombination of the heavier ion, and 
reactions with Hn and Hell to ionisation. 

The electron captured during charge transfer recombination of an oxygen ion (for 
instance Ovn, Ovm) is usually captured in an intermediate quantum state (principal 
quantum number n = 4 — 6). This leads to enhanced line emission from that level 
as compared to the emission from other principal levels, and this signature can be 
used to show the presence of charge transfer reactions. Another signature - actually a 
signature for all recombining plasmas - is of course the enhancement of the forbidden 
line relative to the resonance line in the Ovn triplet (Sect. [5T2.2[l . 

An important example is the charge transfer of highly charged ions from the Solar 
wind with the neutral or weakly ionised Geocorona. Whenever the Sun is more active, 
this process may produce enhanced thermal soft X-ra y emission in ad d ition to the 
steady foreground emission from our own Galaxy. See iBhardwai et al. (2006) for a 
review of X-rays from the Solar System. Care should be taken not to confuse this 
temporary enhanced Geocoronal emission with soft excess emission in a background 
astrophysical source. 



4 Ionisation balance 

In order to calculate the X-ray emission or absorption from a plasma, apart from the 
physical conditions also the ion concentrations must be known. These ion concentra- 
tions can be determined by solving the equations for ionisation balance (or in more 
complicated cases by solving the time-dependent equations). A basic ingredient in these 
equations are the ionisation and recombination rates, that we discussed in the previ- 
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ous section. Here we consider three of the most important cases: collisional ionisation 
equilibrium, non-equilibrium ionisation and photoionisation equilibrium. 



4.1 Collisional Ionisation Equilibrium (CIE) 

The simplest case is a plasma in collisional ionisation equilibrium (CIE). In this case 
one assumes that the plasma is optically thin for its own radiation, and that there is 
no external radiation field that affects the ionisation balance. 

Photo-ionisation and Compton ionisation therefore can be neglected in the case 
of CIE. This means that in general each ionisation leads to one additional free elec- 
tron, because the direct ionisation and excitation-autoionisation processes are most 
efficient for the outermost atomic shells. The relevant ionisation processes are colli- 
sional ionisation and excitation-autoionisation, and the recombination processes are 
radiative recombination and dielectronic recombination. Apart from these processes, 
at low temperatures also charge transfer ionisation and recombination are important. 

We define R z as the total recombination rate of an ion with charge z to charge 
z — 1, and I z as the total ionisation rate for charge z to z + 1. Ionisation equilibrium 
then implies that the net change of ion concentrations n z should be zero: 

z > : n z+1 R z+1 - n z Rz + n z -\I z -i - n z I z = (22) 

and in particular for z = one has 

n%Ri = n Io (23) 

(a neutral atom cannot recombine further and it cannot be created by ionisation). Next 
an arbitrary value for no is chosen, and ([23)) is solved: 

ni=n (Io/Ri). (24) 

This is substituted into (|22p which now can be solved. Using induction, it follows that 

n z+1 =n z {I z /R z+1 ). (25) 

Finally everything is normalised by demanding that 

z 

Z2 n z = "element (26) 

2=0 

where n e i ement is the total density of the element, determined by the total plasma 
density and the chemical abundances. 

Examples of plasmas in CIE are the Solar corona, coronae of stars, the hot intra- 
cluster medium, the diffuse Galactic ridge component. Fig. [7] shows the ion fractions 
as a function of temperature for two important elements. 
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Fig. 7 Ion concentration of oxygen ions (left panel) and iron ions (right panel) as a function 
of temperature in a plasma in Collisional Ionisation Equilibrium (CIE). Ions with completely 
filled shells are indicated with thick lines: the He-like ions O VII and Fe XXV, the Ne-like Fe XVII 
and the Ar-like Fe IX; note that these ions are more prominent than their neighbours. 



4.2 Non-Equilibrium Ionisation (NEI) 

The second case that we discuss is non-equilibrium ionisation (NEI). This situation oc- 
curs when the physical conditions of the source, like the temperature, suddenly change. 
A shock, for example, can lead to an almost instantaneous rise in temperature. How- 
ever, it takes a finite time for the plasma to respond to the temperature change, as 
ionisation balance must be recovered by collisions. Similar to the CIE case we assume 
that photoionisation can be neglected. For each element with nuclear charge Z we 
write: 

^ ~n(Z, t) = A(Z, T(t))n(Z, t) (27) 

where n is a vector of length Z + 1 that contains the ion concentrations, and which is 
normalised according to Eqn. 1261 The transition matrix A is a (Z + 1) x (Z+ 1) matrix 
given by 



A = 



f-I i?i 
J -(h + Ri) R 2 
h 



V 



Rz- 
Iz- 




-{1 2 



\ 





-i + Rz- 



Iz- 



Rz 
-RzJ 



We can write the equation in this form because both ionisations and recombinations 
are caused by collisions of electrons with ions. Therefore we have the uniform scaling 
with n c . In general, the set of equations (|27|l must be solved numerically. The time 
evolution of the plasma can be described in general well by the parameter 



U 



n c dt. 



(28) 
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The integral should be done over a co-moving mass element. Typically, for most ions 
equilibrium is reached for U ~ 10 18 m~ 3 s. We should mention here, however, that the 
final outcome also depends on the temperature history T(t) of the mass element, but 
in most cases the situation is simplified to T(t) — constant. 



4.3 Photoionisation Equilibrium (PIE) 

The third case that we treat are the photoionised plasmas. Usually one assumes equi- 
librium (PIE) , but there are of course also extensions to non-equilibrium photo- ionised 
situations. Apart from the same ionisation and recombination processes that play a 
role for plasmas in NEI and CIE, also photoionisation and Compton ionisation are rel- 
evant. Because of the multiple ionisations caused by Auger processes, the equation for 
the ionisation balance is not as simple as (|22|l . because now one needs to couple more 
ions. Moreover, not all rates scale with the product of electron and ion density, but the 
balance equations also contain terms proportional to the product of ion density times 
photon density. In addition, apart from the equation for the ionisation balance, one 
needs to solve simultaneously an energy balance equation for the electrons. In this en- 
ergy equation not only terms corresponding to ionisation and recombination processes 
play a role, but also several radiation processes (Bremsstrahlung, line radiation) or 
Compton scattering. The equilibrium temperature must be de termined i n an i terative 
way. A classical paper describing such photoionised plasmas is iKallmanl l|l982l ). 



5 Emission processes 

5.1 Continuum emission processes 

In any plasma there are three important continuum emission processes, that we briefly 
mention here: Bremsstrahlung, free-bound emission and two-photon emission. 



5.1.1 Bremsstrahlung 

Bremsstrahlung is caused by a collision between a free electron and an ion. The emis- 
sivity eg (photons m~ 3 s _1 J -1 ) can be written as: 

2V2aa T cn c niZ^ B / m c c 2 \h -E/kT 



es = — 7^ — n^J 5ffe ' (29) 

where a is the fine structure constant, <tt the Thomson cross section, n e and nj the 
electron and ion density, and E the energy of the emitted photon. The factor gg is the 
so-called Gaunt factor and is a dimensionless quantity of order unity. Further, Z e fi is 
the effective charge of the ion, defined as 



where _Ejj is the ionisation energy of hydrogen (13.6 eV), I r the ionisation potential of 
the ion after a recombination, and n r the corresponding principal quantum number. 
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It is also possible to write (|29[) as eg = Pgn e ni with 



3.031 x 10- 21 Z G 2 ff ff off e- g / kr 

^ff = 1/2 ' ( S > 

where in this case Pg is in photons x m s keV - and -BkcV is the energy in keV. The 
total amount of radiation produced by this process is given by 

w 4.856 x lQ- 37 Wm 3 7 2 -E/kT AP 

Wtot = == / Z cB g s e I d£ ke y (32) 

V-'keV J 
n 



From (|29l) we see immediately that the Bremsstrahlung spectrum (expressed in Wm~ 3 
keV -1 ) is flat for E <C kT, and for E > kT it drops exponentially. In order to measure 
the temperature of a hot plasma, one needs to measure near E ~ kT. The Gaunt factor 
gg can be calculated analytically; there are both tables and asymptotic approximations 
available. In general, gg depends on both E/kT and kT/Z cB . 

For a plasma (I29p needs to be summed over all ions that are present in order 
to calculate the total amount of Bremsstrahlung radiation. For cosmic abundances, 
hydrogen and helium usually yield the largest contribution. Frequently, one defines an 
average Gaunt factor Gg by 

Off = { — ) Z ett,i9ef£,i- (33) 



5.1.2 Free-bound emission 

Free-bound emission occurs during radiative recombination (Sect. [3T3. 1 [) . The energy 
of the emitted photon is at least the ionisation energy of the recombined ion (for 
recombination to the ground level) or the ionisation energy that corresponds to the 
excited state (for recombination to higher levels). From the recombination rate (see 
Sect. I3.3TT]) the free-bound emission is determined immediately: 



£fb = ^2 n en\Rr. (34) 
i 

Also here it is possible to define an effective Gaunt factor Gfo. Free-bound emission 
is in practice often very important. For example in CIE for kT = 0.1 keV, free-bound 
emission is the dominant continuum mechanism for E > 0.1 keV; for kT = 1 keV it 
dominates above 3 keV. For kT 3> 1 keV Bremsstrahlung is always the most important 
mechanism, and for kT 0.1 keV free-bound emission dominates. See also Fig. [8] 

Of course, under conditions of photoionisation equilibrium free-bound emission is 
even more important, because there are more recombinations than in the CIE case 
(because T is lower, at comparable ionisation). 
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5.1.3 Two photon emission 

This process is in particular important for hydrogen-like or helium-like ions. After a 
collision with a free electron, an electron from a bound Is shell is excited to the 2s shell. 
The quantum-mechanical selection rules do not allow that the 2s electron decays back 
to the Is orbit by a radiative transition. Usually the ion will then be excited to a higher 
level by another collision, for example from 2s to 2p, and then it can decay radiatively 
back to the ground state (Is). However, if the density is very low (n e <C n e cr it, Ean.l36f- 
I37[) . the probability for a second collision is very small and in that case two-photon 
emission can occur: the electron decays from the 2s orbit to the Is orbit while emitting 
two photons. Energy conservation implies that the total energy of both photons should 
equal the energy difference between the 2s and Is level (-B2phot — ^ls ~ ^2s)- From 
symmetry considerations it is clear that the spectrum must be symmetrical around 
E = 0.5-B 2p hot! and further that it should be zero for E = and E = i^phot- An 
empirical approximation for the shape of the spectrum is given by: 



F(E) ~ yJ S m(irE/E 2phot ). (35) 

An approximation for the critical density below which two photon emission is important 
can be obtained from a co mparison of the radi ative and collisional rates from the upper 
(2s) level, and is given by ( Mewe et al.l . [l98r3 ): 



H - like : n e , crit = 7 x 10 9 m~ 3 Z 9 5 (36) 
He -like: n 0jCrit = 2 x 10 11 m" 3 (Z-l) 95 . (37) 

For example for carbon two photon emission is important for densities below 10 17 m~ 3 , 
which is the case for many astrophysical applications. Also in this case one can deter- 
mine an average Gaunt factor G^phot by averaging over all ions. Two photon emission is 
important in practice for 0.5 < kT < 5 keV, and then in particular for the contributions 
of C, N and O between 0.2 and 0.6 keV. See also Fig. |S] 



5.2 Line emission processes 

Apart from continuum radiation, line radiation plays an important role for thermal 
plasmas. In some cases the flux, integrated over a broad energy band, can be completely 
dominated by line radiation (see Fig. [8j) . The production process of line radiation can 
be broken down in two basic steps: the excitation and the emission process. 

5.2.1 Excitation process 

An atom or ion must first be brought into an excited state before it can emit line 
photons. There are several physical processes that can contribute to this. 

The most important process is usually collisional excitation (Sect. [3~T|) . in particular 
for plasmas in CIE. The collision of an electron with the ion brings it in an excited 
state. 

A second way to excite the ion is by absorbing a photon with the proper energy. 
We discuss this process in more detail in Sect. [6] 

Alternatively, inner shell ionisation (either by the collision with a free electron, 
Sect. l3.2"T1 or by the photoelectric effect, Sect. I3.2T2]| brings the ion in an excited state. 
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* kT=0.04 keV 




Fig. 8 Emission spectra of plasmas with solar abundances. The histogram indicates the total 
spectrum, including line radiation. The spectrum has been binned in order to show better the 
relative importance of line radiation. The thick solid line is the total continuum emission, the 
thin solid line the contribution due to Bremsstrahlung, the dashed line free-bound emission 
and the dotted line two-photon emission. Note the scaling with 

Ee E/kT 

along the y-axis. 



Finally, the ion can also be brought in an excited state by capturing a free elec- 
tron in one of the free energy levels above the ground state (radiative recombination, 
Sect. I3.3TT]) . or through dielectronic recombination ( Sect [533)1 . 



5.2.2 Line emission 



It does not matter by whatever process the ion is brought into an excited state j, 
whenever it is in such a state it may decay back to the ground state or any other lower 
energy level i by emitting a photon. The probability per unit time that this occurs 
is given by the spontaneous transition probability Am (units: s _1 ) which is a number 
that is different for each transition. The total line power Pij (photons per unit time 
and volume) is then given by 

= AijTij (38) 

where nj is the number density of ions in the excited state j. For the most simple case 
of excitation from the ground state g (rate S g j ) followed by spontaneous emission, one 
can simply approximate n g n e S g j = njA g j. From this equation, the relative population 
nj/n g 1 is determined, and then using (|38|l the line flux is determined. In realistic 
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situations, however, things are more complicated. First, the excited state may also 
decay to other intermediate states if present, and also excitations or cascades from 
other levels may play a role. Furthermore, for high densities also collisional excitation 
or de-excitation to and from other levels becomes important. In general, one has to 
solve a set of equations for all energy levels of an ion where all relevant population and 
depopulation processes for that level are taken into account. For the resulting solution 
vector rij, the emitted line power is then simply given by Eqn. (1381) . 

Note that not all possible transitions between the different energy levels are allowed. 
There are strict quantum mechan ical selecti o n rule s that govern which lines are allowed; 
see for instance iHerzberel 1 1944h or iMewel ( 19991 ). Sometimes there are higher order 



processes that still allow a forbidden transition to occur, albeit with much smaller 
transition probabilities A^j . But if the excited state j has no other (fast enough) way 
to decay, these forbidden lines occur and the lines can be quite strong, as their line 
power is essentially governed by the rate at which the ion is brought into its excited 
state j. 

One of the most well known groups of lines is the He- like ls-2p triplet. Usually 
the strongest line is the resonance line, an allowed transition. The forbidden line has 
a similar strength as the resonance line, for the low density conditions in the ISM 
and intracluster medium, but it can be relatively enhanced in recombining plasmas, or 
relatively reduced in high density plasmas like stellar coronal loops. In between both 
lines is the intercombination line. In fact, this intercombination line is a doublet but 
for the lighter elements both components cannot be resolved. But see Fig. HJ] for the 
case of iron. 



5.2.3 Line width 

For most X-ray spectral lines, the line profile of a line with energy E can be approx- 
imated with a Gaussian exp (-AE 2 /2a 2 ) with a given by a/E = av/c where the 
velocity dispersion is 

al = a 2 + kTi/mi. (39) 

Here 7] is the ion temperature (not necessarily the same as the electron temperature), 
and at is the root mean squared turbulent velocity of the emitting medium. For large 
ion temperature, turbulent velocity or high spectral resolution this line width can be 
measured, but in most cases the lines are not resolved for CCD type spectra. 



5.2.4 Resonance scattering 

Resonance scattering is a process where a photon is absorbed by an atom and then 
re-emitted as a line photon of the same energy into a different direction. As for strong 
resonance lines (allowed transitions) the transition probabilities Aij are large, the time 
interval between absorption and emission is extremely short, and that is the reason 
why the process effectively can be regarded as a scattering process. We discuss the 
absorption properties in Sect. 16.31 and have already discussed the spontaneous emission 
in Sect. [5^21 

Resonance scattering of X-ray photons is potentially important in the dense cores 
of some clusters of galaxies for a few transitions of the most abundant elements, as 
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first shown bv lGn'fanov et all ( 19871 '). The optical depth for scattering can be written 
conveniently as (cf. also Sect. 16.3]) : 



4240 /iV 2 4 



( nj_\ (nz\ ( M \ 



r 0.0522M V ? nn 1 1/2 

^kcV 1 1 + 



1/2 

— , (40) 



TkeV 

where / is the absorption oscillator strength of the line (of order unity for the strongest 
spectral lines), -EkeV the energy in keV, N24 the hydrogen column density in units of 
10 24 m -2 , rij the number density of the ion, ni the number density of the element, 
M the atomic weight of the ion, T^y the ion temperature in keV (assumed to be 
equal to the electron temperature) and D100 the micro-turbulence velocity in units 
of 100 km/s. Resonance scattering in clusters causes the radial intensity profile on 
the sky of an emission line to become weaker in the cluster core and stronger in the 
outskirts, without destroying photons. By comparing the radial line profiles of lines 
with different optical depth, for instance the Is— 2p and Is— 3p lines of O VII or Fexxv, 
one can determine the optical depth and hence constrain the amount of turbulence in 

the cluster. 

Another important application was proposed bv lChurazov et all (|200lh . They show 
that for WHIM filaments the resonance line of the O VII triplet can be enhanced signif- 
icantly compared to the thermal emission of the filament due to resonant scattering of 
X-ray background photons on the filament. The ratio of this resonant line to the other 
lines of the triplet therefore can be used to estimate the column density of a filament. 



5.3 Some important line transitions 

In Tables [4j{5] we list the 100 strongest emission lines under CIE conditions. Note that 
each line has its peak emissivity at a different temperature. In particular some of the 
H-like and He-like transitions are strong, and further the so-called Fe-L complex (lines 
from n — 2 in Li-like to Ne-like iron ions) is prominent. At longer wavelengths, the 
L-complex of Ne Mg, Si and S gives strong soft X-ray lines. At very short wavelengths, 
there are not many strong emission lines: between 6-7 keV, the Fe-K emission lines are 
the strongest spectral features. 



6 Absorption processes 

6.1 Continuum versus line absorption 

X-rays emitted by cosmic sources do not travel unattenuated to a distant observer. This 
is because intervening matter in the line of sight absorbs a part of the X-rays. With 
low-resolution instruments, absorption can be studied only through the measurement 
of broad-band flux depressions caused by continuum absorption. However, at high 
spectral resolution also absorption lines can be studied, and in fact absorption lines 
offer more sensitive tools to detect weak intervening absorption systems. We illustrate 
this in Fig. [9] At an Ovni column density of 10 21 m -2 , the absorption edge has an 
optical depth of 1 %; for the same column density, the core of the Lya line is already 
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Table 4 The strongest emission lines for a plasma with proto-solar abundances (jLoddcrs, 
2003) in the X-ray band 43 A < A < 100 A. At longer wavelengths sometimes a few lines from 
the same multiplct have been added. All lines include unresolved dielectronic satellites. T max 
(K) is the temperature at which the emissivity peaks, Q ma x = P/(n e nH)i with P the power 
per unit volume at T max , and Q ma x is in units of 10 — 36 W m 3 . 
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Table 5 As Table H but for A < 43 A. 
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saturated and even for 100 times lower column density, the line core still has an optical 
depth of 5 %. 




6.2 Continuum absorption 

Continuum absorption can be calculated simply from the photoionisation cross sections, 
that we discussed already in Sect. 13.2.21 The total continuum opacity T con t can be 
written as 

r cont (E) = TVhCW (E) = NiCTi (E) (41) 

i 

i.e, by averaging over the various ions i with column density TV,;. Accordingly, the 
continuum transmission T(E) of such a clump of matter can be written as T(E) = 
cxp (— r con t (E)). For a worked out example see also Sect. 16.61 



6.3 Line absorption 

When light from a background source shines through a clump of matter, a part of 
the radiation can be absorbed. We discussed already the continuum absorption. The 
transmission in a spectral line at wavelength A is given by 

T(A) = e~ r ( A ) (42) 

with 

r(A) = to<^(A) (43) 
where <fi{\) is the line profile and tq is the opacity at the line centre Aq, given by: 
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Apart from the fine structure constant a and Planck's constant h, the optical depth 
also depends on the properties of the absorber, namely the ionic column density TV, 
and the velocity dispersion ov. Furthermore, it depends on the oscillator strength / 
which is a dimensionless quantity that is different for each transition and is of order 
unity for the strongest transitions. 

In the simplest approximation, a Gaussian profile ip(X) = exp — (A — Xo) 2 /b 2 ) can 
be adopted, corresponding to pure Doppler broadening for a thermal plasma. Here 
b = \f2<7 with a the normal Gaussian root-mean-square width. The full width at half 
maximum of this profile is given by Vm256cr or approximately 2.35o. One may even 
include a turbulent velocity at into the velocity width ov, such that 

a 2 = at + kT/m, (45) 

with m.j the mass of the ion (we have tacitly changed here from wavelength to velocity 
units through the scaling Z\A/Ao = Av/c). 

The equivalent width W of the line is calculated from 

oo 

W= ^~ I [ 1 - ex P(- r o e_!/2/2 )] d y- ( 46 ) 

— oo 

For the simple case that to <C 1, the integral can be evaluated as 

T0 « 1 : W= ~^tf = 2 a ^ fX Ni - (47) 

A Gaussian line profile is only a good approximation when the Doppler width of 
the line is larger than the natural width of the line. The natural line profile for an 
absorption line is a Lorentz profile tp(X) = 1/(1 + x 2 ) with x — 4-kAv/A. Here Av is 
the frequency difference v — v$ and A is the total transition probability from the upper 
level downwards to any level, including all radiative and Auger decay channels. 

Convolving the intrinsic Lorentz profile with the Gaussian Doppler profile due to 
the thermal or turbulent motion of the plasma, gives the well-known Voigt line profile 

<p = H(a,y) (48) 

where 

a = AX /4-Kb (49) 

and y = cA\/b\. The dimensionless parameter a (not to be confused with the a in 
Eqn.[4]| represents the relative importance of the Lorentzian term (~ A) to the Gaussian 
term (~ b). It should be noted here that formally for a > the parameter tq defined 
by (144 1) is not exactly the optical depth at line centre, but as long as a <C 1 it is a fair 
approximation. For a ^> 1 it can be shown that H(a,0) — > l/ay/w. 

6.4 Some important X-ray absorption lines 

There are several types of absorption lines. The most well-known are the normal strong 
resonance lines, which involve electrons from the outermost atomic shells. Examples 
are the ls-2p line of Ovu at 21.60 A, the Ovm Lya doublet at 18.97 A and the well 
known 2s-2p doublet of Ovi at 1032 and 1038 A. See lRichter et alJboOSl - Chapter 3, 
this volume, for an extensive discussion on these absorption lines in the WHIM. 



27 



The other class are the inner-shell absorption lines. In this case the excited atom is 
often in an unstable state, with a large probability for the Auger process (Sect. l372,4|l . 
As a result, the parameter a entering the Voigt profile (Eqn. |49)) is large and therefore 
the lines can have strong damping wings. 

In some cases the lines are isolated and well resolved, like the O VII and O VIII ls- 
2p lines mentioned above. However, in other cases the lines may be strongly blended, 
like for the higher principal quantum number transitions of any ion in case of high 
column densities. Another important group of strongly blended lines are the inner- 
shell transitions of heavier elements like iron. They form so-called unresolved transition 
arrays (UTAs); the individual lines are no longer rec ognisable. The fir st detection of 
these UTAs in astrophysical sources was in a quasar ijSako et all l200lh . 

In Table[6]we list the 70 most important X-ray absorption lines for A < 100 A. The 
importance was determined by calculating the maximum ratio W/X that these lines 
reach for any given temperature. 

Note the dominance of lines from all ionisation stages of oxygen in the 17-24 A 
band. Furthermore, the Feix and Fexvn lines are the strongest iron features. These 
lines are weaker than the oxygen lines because of the lower abundance of iron; they 
are stronger than their neighbouring iron ions because of somewhat higher oscillator 
strengths and a relatively higher ion fraction (cf. Fig. [7} . 

Note that the strength of these lines can depend on several other parameters, like 
turbulent broadening, deviations from CIE, saturation for higher column densities, etc., 
so some care should be taken when the strongest line for other physical conditions are 
sought. 



6.5 Curve of growth 

In most practical situations, X-ray absorption lines are unresolved or poorly resolved. 
As a consequence, the physical information about the source can only be retrieved from 
basic parameters such as the line centroid (for Doppler shifts) and equivalent width W 
(as a measure of the line strength). 

As outlined in Sect. 16.31 the equivalent width W of a spectral line from an ion i 
depends, apart from atomic parameters, only on the column density iVj and the velocity 
broadening a v . A useful diagram is the so-called curve of growth. Here we present it 
in the form of a curve giving W versus iVj for a fixed value of a v . 

The curve of growth has three regimes. For low column densities, the optical depth 
at line centre to is small, and W increases linearly with iVj, independent of a v (Eqn.HT)). 
For increasing column density, tq starts to exceed unity and then W increases only log- 
arithmically with ATj. Where this happens exactly, depends on <r v . In this intermediate 
logarithmic branch W depends both on a v and iVj, and hence by comparing the mea- 
sured W with other lines from the same ion, both parameters can be estimated. At 
even higher column densities, the effective opacity increases faster (proportional to 
\/77i), because of the Lorentzian line wings. Where this happens depends on the Voigt 
parameter a (Eqn. [49}. In this range of densities, W does not depend any more on a v . 

In Fig. [10] we show a few characteristic examples. The examples of O I and O VI 
illustrate the higher W for UV lines as compared to X-ray lines (cf. Eqn. I47[l as well 
as the fact that inner shell transitions (in this case the X-ray lines) have larger a- 
values and hence reach sooner the square root branch. The example of O VII shows 
how different lines from the same ion have different equivalent width ratios depending 
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Table 6 The most important X-ray absorption lines with A < 100 A. Calculations are done for 
a plasma in CIE with proto-solar abundances (Loddcrs. 2003) and only thermal broadening. 
The calculations are done for a hydrogen column density of 10 24 m -2 , and we list equivalent 
widths for the temperature T max (in K) where it reaches a maximum. Equivalent widths are 
calculated for the listed line only, not taking into account blending by other lines. For saturated, 
overlapping lines, like the Ovm Lya doublet near 18.97 A, the combined equivalent width can 
be smaller than the sum of the individual equivalent widths. 
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on ov, hence offer an opportunity to determine that parameter. The last frame shows 
the (non)similarties for similar transitions in different ions. The innershell nature of 
the transition in O vi yields a higher a value and hence an earlier onset of the square 
root branch. 
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Fig. 10 Equivalent width versus column density for a few selected oxygen absorption lines. 
The curves have been calculated for Gaussian velocity dispersions a = b/\/2 of 10, 50 and 
250 kms" 1 , from bottom to top for each line. Different spectral lines are indicated with 
different line styles. 



6.6 Galactic foreground absorption 

All radiation from X-ray sources external to our own Galaxy has to pass through the 
interstellar medium of our Galaxy, and the intensity is reduced by a factor of e~ T ^ E ^ 
with the optical depth t(E) = Ci(E) J n,i(£)d£, with the summation over all relevant 
ions i and the integration over the line of sight d£. The absorption cross section <Ti{E) 
is often taken to be simply the (continuum) photoionisation cross section, but for 
high-resolution spectra it is necessary to include also the line opacity, and for very 
large column densities also other processes such as Compton ionisation or Thomson 
scattering. 

For a cool, neutral plasma with cosmic abundances one often writes 

r = a cB (E)N H (50) 

where the hydrogen column density TVh = / n^dx. In a e ff(E) all contributions to the 
absorption from all elements as well as their abundances are taken into account. 

The relative contribution of the elements is also made clear in Fig. 1111 Below 
0.28 keV (the carbon edge) hydrogen and helium dominate, while above 0.5 keV in 
particular oxygen is important. At the highest energies, above 7.1 keV, iron is the 
main opacity source. 
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Fig. 11 Left panel: Neutral interstellar absorption cross section per hydrogen atom, scaled 
with . The most important edges with associated absorption lines are indicated, as well as 
the onset of Thompson scattering around 30 keV. Right panel: Contribution of the various 
elements to the total absorption cross section as a function of energy, for solar abundances. 




0.01 0.1 1 10 100 0.01 0.1 1 10 100 

Energy (keV) Energy (keV) 

Fig. 12 Left panel: Column density for which the optical depth becomes unity. Right panel: 
Distance for which r = 1 for three characteristic densities. 



Yet another way to represent these data is to plot the column density or distance for 
which the optical depth becomes unity (Fig. [12} • This figure shows that in particular 
at the lowest energies X-rays are most strongly absorbed. The visibility in that region 
is thus very limited. Below 0.2 keV it is extremely hard to look outside our Milky Way. 

The interstellar medium is by no means a homogeneous, neutral, atomic gas. In 
fact, it is a collection of regions all with different physical state and composition. This 
affects its X-ray opacity. We briefly discuss here some of the most important features. 

The ISM contains cold gas (< 50 K), warm, neutral or lowly ionised gas (6 000 — 
10 000 K) as well as hotter gas (a few million K). We have seen before (see Sect. I3.2T2"]! 
that for ions the absorption edges shift to higher energies for higher ionisation. Thus, 
with instruments of sufficient spectral resolution, the degree of ionisation of the ISM 
can be deduced from the relative intensities of the absorption edges. Cases with such 
high column densities are rare, however (and only occur in some AGN outflows) , and for 
the bulk of the ISM the column density of the ionised ISM is low enough that only the 
narrow absorption lines are visible (see Fig. [9}. These lines are only visible when high 
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o 4U 1820-303 




Fig. 13 Left panel: Simulated 100 ks absorption spectrum as observed with the Explorer of 
Diffuse Emission and Gamma-ray Burst Explosions (EDGE), a mission proposed for ESA's 
Cosmic Vision program. Th e para meters of the simulated source are similar to those of 
4U 1820 - 303 l lYao fc Wanel . 120061) . The plot shows the residuals of the simulated spectrum 
if the absorption lines in the model are ignored. Several characteristic absorption features of 
both neutral and ionised gas are indicated. 

Right panel: Simulated spectrum for the X-ray binary 4U 1820 - 303 for 100 ks with the WFS 
instrument of EDGE. The simulation was done for all absorbing oxygen in its pure atomic state, 
the models plotted with different line styles show cases where half of the oxygen is bound in 
CO, water ice or olivine. Note the effective shift of the absorption edge and the different fine 
structure near the edge. All of this is well resolved by EDGE, allowing a determination of the 
molecular composition of dust in the line of sight towards this source. The absorption line at 
0.574 is due to highly ionised Ovn. 



spectral resolution is used (see Fig. I13[l . It is important to recognise these lines, as they 
should not be confused with absorption line from within the source itself. Fortunately, 
with high spectral resolution the cosmologically redshifted absorption lines from the 
X-ray source are separated well from the foreground hot ISM absorption lines. Only for 
lines from the Local Group this is not possible, and the situation is here complicated 
as the expected temperature range for the diffuse gas within the Local Group is similar 
to the temperature of the hot ISM. 

Another ISM component is dust. A significant fracti on of some atoms ca n be bound 
in dust grains with varying sizes, as shown below (from lwilms et al"1l2000l ): 

H He C N O Ne Mg Si S Ar Ca Fe Ni 
0.5 0.4 0.8 0.9 0.4 0.997 0.7 0.96 

The numbers represent the fraction of the atoms that are bound in dust grains. Noble 
gases like Ne and Ar are chemically inert hence are generally not bound in dust grains, 
but other elements like Ca exist predominantly in dust grains. Dust has significantly 
different X-ray properties compared to gas or hot plasma. First, due to the chemical 
binding, energy levels are broadened significantly or even absent. For example, for oxy- 
gen in most bound forms (like H2O) the remaining two vacancies in the 2p shell are 
effectively filled by the two electrons from the other bound atom(s) in the molecule. 
Therefore, the strong ls-2p absorption line at 23.51 A (527 eV) is not allowed when 
the oxygen is bound in dust or molecules, because there is no vacancy in the 2p shell. 
Transitions to higher shells such as the 3p shell are possible, however, but these are 
blurred significantly and often shifted due to the interactions in the molecule. Each 
constituent has its own fine structure near the K-edge (Fig. 113b ). This fine structure 
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offers therefore the opportunity to study the (true) chemical composition of the dust, 
but it should be said that the details of the edges in different important compounds are 
not always (accurately) known, and can differ depending on the state: for example wa- 
ter, crystalline and amorphous ice all have different characteristics. On the other hand, 
the large scale edge structure, in particular when observed at low spectral resolution, 
is not so much affected. For sufficiently high column densities of dust, self-shielding 
within the grains should be taken into account, and this reduces the average opacity 
per atom. 

Finally, we mention here that dust also causes scattering of X-rays. This is in 
particular important for higher column densities. For example, for the Crab nebula 
(N H = 3.2 x 10 25 m~ 2 ), at an energy of 1 keV about 10 % of all photons are scattered 
in a halo, of which the widest tails have been observed out to a radius of at least half 
a degree; the scattered fraction increases with increasing wavelength. 



7 Galactic foreground emission 

The interstellar medium of our Galaxy is a complex medium. While all phases of 
the ISM can be seen in X-ray absorption towards bright X-ray sources (see previous 
section), only the hot phase is seen in emission and contributes to the cosmic X-ray 

background. 

This cosmic X-ray background has different components. iKuntz fc Snowdenl ( 200 lh 



distinguish four different components as outlined below. First, there is an absorbed 
power-law like component, consisting mostly of unresolved extragalactic point sources. 
With high spatial resolution like for example available on Chandra, a large part of this 
component can be resolved into the individual point sources. The second component is 
the so-called Local Hot Bubble, most likely a large, old supernova remnant embedding 
our Solar system. The million degree gas of this diffuse medium is almost unabsorbed, 
as it is nearby and there is little intervening neutral gas. Finally, there are the soft and 
hard diffuse components, consisting of absorbed, diffuse thermal emission which arises 
from the disk of our Galaxy but may also contain contributions from a Galactic halo 
and distant WHIM emission. 

Fig. 1141 shows a simulated spectrum of the total X-ray background. Note however 
that the details of this spectrum may differ significantly for different parts of the sky. 
There are variations on large scales, from the Galactic plane to the poles. There are 
also smaller scale variations. For instance, there are large loops corresponding to old 
supernova remnants or (super)bubbles, and large molecular cloud complexes obscuring 
the more distant background components. 



8 Cooling function 

The cooling of a plasma by thermal radiation is usually expressed by the cooling 
function A(T). This is essentially the bolometric power P emitted as thermal radiation 
by the plasma, normalised by the emission measure Y = n c nnV . Accordingly, the 
characteristic cooling time, defined as the thermal energy of the plasma divided by the 
emitted power, t coo i = |(1 + n,[/n e )kT/ Ann. 

Cooling is important for the evolution of several plasmas such as the cool cores of 
clusters but also for the global evolution of the WHIM. A set of useful cooling curves are 
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Fig. 14 Typical Cosmic X-ray background spectrum. The thick solid line indicates the total 
spectrum, the dashed line the absorbed power law contribution due to unresolved point sources. 
The contribution of the local hot bubble is indicated by the thin, grey line. It dominates the 
line emission below ~0.2 keV. The other thermal components are not indicated separately, 
but they are included in the total; they contribute most of the line flux above 0.2 keV. The 
calculation is done with a spectral resolution of 1 eV. 
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Fig. 15 Cooling curves A(T) for plasmas of different composition, under CIE conditions. From 
top to bottom: 1, 0.3, 0.1, 0.01 and 0.001 times solar abundances. After [Suthcr land & Dopita 
<|1993T >. 



given bv lSutherland fc Dopital (|l993h . An example is shown in Fig. 1151 These are valid 
under CIE conditions. When photoionisation is important, modifications are needed. 
Here we have normalised A to the emission measure of the plasma, but note that 
sometimes other scalings are being used (for example expressions containing the ion 
d ensity n\ or rip). Note that the Solar abundances used in Fig. [15] are the older values 
of lAnders fc Grevessel (|l98Sh which differ significantly from more recent estimates, for 
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example iLodderj (|2003h . in particular for important coolants such as iron (new values 
26 % lower) and oxygen (new values 32 % lower). The big hump at low temperatures 
is produced mainly by line radiation, while the high temperature ~ j 10 - 5 tail above 
10 7 K is mainly produced by Bremsstrahlung. 

Finally, a more recent set of cooling cur ves, suited in parti cular for lower temper- 
atures and primordial gas, was published bv lMaio et al. | 120071 ). 



9 Ionisation and excitation by non-thermal electrons 

In most cases, astrophysical plasmas are modelled using Maxwellian velocity distribu- 
tions, even under NEI or PIE conditions. Often this is justified, but there are situations 
where deviations from a Maxwellian distribution occur. Examples are accelerated elec- 
trons penetrating cooler regions, as may occur in or near shocks or flares. We refrain 
here from a full discussion, but only mention here that spectral codes such as SPEX 
(Sect. [Ten")) allow the calculation of the ionisation balance and the emanating spectrum 
under such circumstances, provided the true electron distribution can be approximated 
by the sum of non-relativistic Maxwellians. When there are strong relativistic electron 
contributions, this approach cannot be used anymore as here relativistic corrections to 
the cross sections are needed, which are usually not taken into account in most plasma 
codes. 



10 Plasma modelling 

10.1 Available codes 

There are several spectral codes available that can be used to model hot plasmas. 
Without being com plete, we mentio n here briefly the following codes. For a more 
extensive review, see iRavmondl l|2005l l. 



First, there are codes which calculate predominantly spectra in Collisional Ionisa- 
tion Equilibrium: 

1. CHIANT^, (named after the wine region near F i renze , lLandi et alJliood ) evolved 
from early work by lLandini fc Monsignori Fossil (|l970h : originally, it focussed on 



the calculation of Solar spectra and hence EUV lines are well represented in this 

code. _ _ 

APEC0 (As trophysical Plasma Emiss ion Code, 1 Smith et all l200l] ) evolved from 



the original IRavmond fc Smith] (jl987t i code; originally, the main goal was to cal- 



culate the emissivity and cooling curve properly . This model is available in the 



culate tne emissivity and cooling curve pr operly, 
popular fitting package XSPEC0, (|Arnaudl. Il99rji. 



Mekal (after MEwe, KAastra & Liedahl; see Mewe et al.l |l99ll) evolved from early 
work by Mewe and collaborators starting from 1970 onwards. The Mekal model is 
also incorporated in XSPEC. 



http:/ /www. arcetri.astro.it/science/chianti/chianti. html 

2 http://cxc.harvard.edu/atomdb 

3 http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspec 
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The Mekal model is now part of the SPEX package^, ijKaastra et allll996bh . SPEX 
is a spectral fitting package developed at SRON Utrecht based on the original CIE 
models of Mewe and collaborators. It has grown significantly and also contains NEI 
and a few simple PIE models; more models are under development. The package offers 
the opportunity to plot or tabulate also various physical parameters of interest of the 
model, apart from the spectrum, and the atomic physics used in the Mekal model 
has been significantly updated in the SPEX version as compared to the older version 
available in XSPEC. A majority of all plots shown in this paper were produced using 
SPEX. 

There are also models focussing upon photo-ionised plasmas. We mention here two 
codes (apart from SPEX): 

1. XSTAE0 evolved from the original work of iKallmanl (|l982l ). focussing upon the 
X-ray band. It is relatively straightforward to use and there are options to use it 
with XSPEC. 

2. Cloud j0 ([Ferland et alJ . ll998h evolved from work by Ferland and collaborators that 
originally was meant to model optical emission line regions. It gradually evolved 
to the X-ray band. It is somewhat more complicated to use, but has more options 
than XSTAR. 



10.2 Important parameters 

In order to calculate X-ray spectra for collisional plasmas, there is a set of parameters 
that determine the spectral shape and flux. The two most important parameters are 
the electron temperature T c and emission measure Y of the gas. The emission measure 
is sometimes not well or even incorrectly defined. We follow here the definition used 
in SPEX, namely Y = f n e nndV where n e and njj are the electron and hydrogen 
density (whether ionised or not). Sometimes people use "n 2 " in this expression, but 
it is not always obvious whether ion density, baryon density, electron density, total 
particle density or something else is meant. The next important parameters are the 
abundances of the elements; these do not only determine the line emission spectrum, 
but also affect the continuum emission (Sect. 15.1 p . Note that the abundances also 
affect the precise values of rte/n n and ni/nu- For a fully ionised plasma with proto- 
solar abundances jLoddersL boost ) these ratios are 1. 201 and 1.097, respectively (for the 
older lAnders fc Grevessel (|l989h abundances these ratios are 1.209 and 1.099). 

For ionising plasmas (NEI) also the parameter U = J n e dt (see Eqn. I28|l is impor- 
tant, as this describes the evolution of a shocked or recombining gas towards equilib- 
rium. 

When high spectral resolution observations are available, other parameters are 
important, like the ion temperature and turbulent velocity. These parameters deter- 
mine the width of the spectral lines (see Sect. [572721 and Sect. [673}. Note that in non- 
equilibrium plasmas, the electron and ion tempera tures are no t neces sarily the same. In 
an extreme case, the supernova remnant SN 1006, IVink et all (|2003h measured an oxy- 



gen ion temperature that is 350 times larger than the electron temperature of 1.5 keV. 
SPEX has also options to approximate the spectra when there are non-Maxwellian tails 



4 http://www.sron.nl/spex 

5 http: / /heasarc. gsfc.nasa.gov/docs/software/xstar/xstar. html 

6 http://www.nublado.org 
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to the electron distribu tion (Sect. El). These t ails may be produced by non-thermal ac- 
celeration mechanisms i Petrosian et al.ll2008l - Chapter 9, this volume). 

In general, for the low densities encountered in the diffuse ISM of galaxies, in 
clusters and in the WHIM, the absolute electron density is not important (except for 
the overall spectral normalisation through the emission measure Y). Only in the higher 
density regions of for example stars X-ray line ratios are affected by density effects, 
but this is outside the scope of the present paper. 

However, if the density becomes extremely low, such as in the lowest density parts 
of the WHIM, the gas density matters, as in such situations photoionisation by the 
diffuse radiation field of galaxies or the cosmic background becomes important. But 
even in such cases, it is not the spectral shape but the ionisation (im)balance that 
is affected most. In those cases essentially the ratio of incoming radiation flux to gas 
density determines the photoionisation equilibrium. 



10.3 Multi- Temperature emission and absorption 



Several astrophysical sources cannot be characterised by a single temperature. There 
may be various reasons for this. It may be that the source has regions with different 
temperatures that cannot be spatially resolved due to the large distance to the source, 
the low angular resolution of the X-ray instrument, or superposition of spectra along 
the line of sight. Examples of these three conditions are stellar coronae (almost point 
sources as seen from Earth) , integrated cluster of galaxies spectra, or spectra along the 
line of sight through the cooling core of a cluster, respectively. But al so in other less 
obvious cases multi-temperature structure can be present. For example iKaastral (2004) 
showed using deprojected spectra that at each position in several cool core clusters there 
is an intrinsic temperature distribution. This may be due to unresolved smaller regions 
in pressure equilibrium with different densities and hence different temperatures. 




Fig. 16 Reconstructed DEM distribution for a simulated spectrum containing two equally 
strong isothermal components of 0.6 and 2.4 keV (left panel) or a uniform distribution between 
0.6 and 2.4 keV (right panel). Thick solid line: the input model; thin solid line: regularisation 
method; dot-dashed line: polynomial method; dashed line: genetic algorithm; dotted line: two 
broadened Gaussian components. The corresponding best-fit spectra differ in \ 2 by less than 
2 for each case. Obviously, small scale detail cannot be recovered but bulk properties such as 
the difference between a bimo dal DEM distribution (left panel) or a flat distribution (right 
panel) can be recovered. From jKaastra et al l ll996al l. See that paper for more details. 
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The most simple solution is to add as many temperature components as are needed 
to describe the observed spectrum, but there are potential dangers to this. The main 
reason is temperature sensitivity. Most spectral temperature indicators (line fluxes, 
Bremsstrahlung exponential cut-off) change significantly when the temperature is changed 
by a factor of two. However, for smaller temperature changes, the corresponding spec- 
tral changes are also minor. For example, a source with two components of equal emis- 
sion measure Y and kT=0.99 and 1.01 keV, respectively, is almost indistinguishable 
from a source with the emission measure 2Y and temperature 1.00 keV. Obviously, 
the most general case is a continuous differential emission measure distribution (DEM) 
D(T) = dY(T)/dT. Typically, in each logarithmic temperature range with a width of 
a factor of two, one can redistribute the differential emission measure D(T) in such 
a way that the total emission measure as integrated over that interval as well as the 
emission-measure weighted temperature remain constant. The resulting spectra are 
almost indistinguishable. This point is illustrated in Fig. 1161 

Th ere are several methods to reconstruct D(T). We refer here to iKaastra et al] 
( 1996al l for more details. All these methods have been implemented in SPEX. They 
have been applied to stell ar coronae or S olar spectra, but recently also to clusters of 
galaxies. See for example iKaastral ( 20041 ) . In the last case, an analytical power-law 
approximation D{T) ~ T 1 / 01 nas been used. 

Also for photo-ionised plasmas a similar technique can be applied. Instead of the 
temperature T, one uses the ionisation parameter £ (XSTAR users) or U (Cloudy 
users). In some cases (like Seyfert 2 galaxies with emission spectra) one recovers D(£) 
from the observed emission spectra; in other cases (like Seyfert 1 galaxies with absorp- 
tion spectra), a continuum model is applied to the observed spectrum, and D(£) is 
recovered from the absorption features in the spectrum. Most of what has been said 
before on D(T) also applies to this case, i.e. the resolution in £ and the various methods 
to reconstruct the continuous distribution D(£). 
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